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RADAR  WAVEFORM  SYNTHESIS  FOR  TARGET  IDENTIFICATION 


Abstract 

A  new  scheme  for  radar  detection  and  discrimination,  the  ra¬ 
dar  waveform  synthesis  method,  is  investigated.  This  scheme  con¬ 
sists  of  synthesizing  an  aspect- independent  waveform  for  the 
incident  radar  signal  which  excites  an  arbitrarily  oriented  tar¬ 
get  in  such  a  way  that  the  return  radar  signal  from  the  target 
contains  only  a  single  natural  resonance  mode  of  the  target  in  the 
late-time  period.  When  the  synthesized  incident  radar  signal  for 
exciting  a  particular  natural  mode  of  a  known,  preselected  target 
is  applied  to  a  wrong  target,  the  return  radar  signal  will  be 
significantly  different  from  that  of  the  expected  natural  mode, 
thus,  the  wrong  target  can  be  sensitively  discriminated.  The 
selection  of  an  optimum  pulse  duration  of  the  required  incident 
signal  for  the  purpose  of  shaping  its  waveform  is  also  studied. 

Three  kinds  of  targets,  an  arbitrarily  oriented  thin  wire,  a 
sphere  and  a  normally  oriented  infinite  cylinder,  have  been  in¬ 
vestigated.  The  cases  of  an  arbitrarily  oriented  wire  and  a  sphere 
have  been  completed  and  the  results  are  reported  in  Parts  1  and  2. 
The  case  of  infinite  cylinder  will  be  completed  shortly  and  in¬ 
complete  results  of  this  case  are  outlined  in  Part  3. 
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1.  Introduction 

In  recent  years,  research  on  radar  target  identification  and 
discrimination  utilizing  a  short-pulse  waveform  has  been  conducted 
by  a  number  of  workers  [1-8].  One  interesting  scheme  is  to  irradi¬ 
ate  a  target  with  a  simple  waveform  such  as  an  impulse,  a  step  or 
a  ramp  signal,  and  then  analyze  the  scattered  field  from  the  tar¬ 
get  in  terms  of  natural  resonance  modes  of  the  target.  It  is  known 
that  the  waveform  of  the  scattered  field  is  aspect  dependent,  but 
the  set  of  natural  resonant  frequencies  extracted  from  the  scattered 
field  is  independent  of  the  aspect  angle  [9-13].  Using  this  pro¬ 
perty,  a  target  can  be  identified  if  the  extracted  set  of  natural 
frequencies  is  compared  with  the  collection  of  known  data  on  the 
natural  frequencies  of  various  targets.  Two  different  targets  can 
also  be  differentiated  if  the  two  sets  of  natural  frequencies  are 
compared.  One  of  the  problems  associated  with  this  scheme  is 
difficulty  in  obtaining  accurate  natural  frequencies  of  the  target 
from  a  noisy  scattered  field. 

In  this  paper,  an  inverse  scheme,  to  be  called  the  radar-wave¬ 
form-synthesis  method,  is  investigated.  Instead  of  analyzing  the 
field  scattered  by  the  target  in  terms  of  its  natural  resonance 
modes,  this  new  scheme  synthesizes  the  waveform  of  the  incident 
radar  signal  in  such  a  way  that,  when  it  excites  the  target,  the 
return  radar  signal  contains  only  a  single  natural  mode  of  the 
target.  It  will  be  shown  that  when  the  incident  radar  signal 
synthesized  to  excite  a  particular  natural  mode  of  a  preselected 
target  is  applied  to  a  different  target,  the  return  signal  will 
be  significantly  different  from  that  of  the  expected  natural  mode. 
The  wrong  target  can  thus  be  sensitively  discriminated. 

The  simplest  case  of  this  radar  waveform  synthesis  scheme  has 
been  studied  by  Chen  [14]  for  the  case  of  a  thin  wire  irradiated 
by  a  radar  pulse  at  normal  incidence.  For  this  case,  it  is  possi¬ 
ble  to  synthesize  a  required  waveform  for  the  incident  radar  wave 
to  excite  a  single-mode  return  response  at  all  post- incidence  times. 
When  this  study  is  generalized  to  oblique  incidence,  difficulties 
are  encountered  in  obtaining  a  realizable  required  incident  wave¬ 
form  for  exciting  a  single-mode,  scattered  field.  Furthermore, 
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the  required  incident  radar  signal  appears  to  be  aspect  dependent. 
This  difficulty  arises  because  there  exists  a  finite  transit  time 
for  an  obliquely-oriented  wire,  i.e.  a  finite  time  for  an  impulse 
to  pass  the  wire.  The  impulse  response  of  this  wire  consists  of 
an  early-time,  forced  response  in  addition  to  the  sum  of  natural 
modes  which  describes  a  normally  oriented  wire.  This  early-time, 
forced  impulse  response  is  difficult  to  approximate  analytically, 
and  consequently  is  responsible  for  problems  encountered  when 
synthesizing  an  incident  radar  signal  to  excite  a  single-mode, 
scattered  field  at  all  post-incidence  times. 

To  overcome  this  difficulty,  we  have  concentrated  on  the  be¬ 
havior  of  the  late-time  response  of  the  wire,  and  have  found  a 
scheme  to  synthesize  the  required  waveform  for  an  incident  radar 
signal  of  finite  duration  to  excite  a  single-mode,  scattered  field 
in  the  late-time  period  (where  the  early-time  impulse  response  is 
not  required,  since  that  period  has  elapsed).  More  significantly, 
this  synthesized  incident  radar  signal  was  found  to  be  aspect  in¬ 
dependent  . 

Initially  the  impulse  response  of  an  obliquely  oriented  wire 
is  approximated  analytically.  This  response  is  found  to  consist 
of  a  forced  response,  which  exists  only  during  the  early-time 
period,  augmenting  a  causal  stun  of  natural  resonance  modes.  It  is 
next  demonstrated  that  an  aspect- independent ,  incident  radar  signal 
which  excites  a  single-mode  return  signal  after  the  early-time 
period  can  be  synthesized.  The  aspect-independent,  required  wave¬ 
forms  of  the  incident  radar  signals  which  excite  various  single¬ 
mode  return  signals  from  an  arbitrarily  oriented  wire  are  then  ob¬ 
tained.  Numerical  examples  on  required  waveforms  for  exciting 
various  single-mode  return  signals,  and  those  on  the  return  signals 
from  a  wire  oriented  at  various  angles  are  given.  Numerical  ex¬ 
amples  are  also  given  to  show  that  when  the  synthesized  incident 
signal  which  excites  a  particular  natural  mode  of  a  preselected 
target  is  applied  to  a  different  target,  the  return  signal  from  the 
wrong  target  becomes  significantly  different  from  that  of  the  ex¬ 
pected  natural  mode.  The  possibility  of  shaping  the  waveform  of 
the  required  incident  signal  by  adjusting  the  duration  of  the 
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incident  signal,  and  other  features  of  the  radar  waveform  synthe¬ 
sis  method  are  also  discussed. 

2.  Geometry  of  Problem 

A  thin-wire  cylinder  of  length  L  and  radius  a  is  illuminated 
by  an  incident-wave  radar  signal  at  angle  e  as  indicated  in  Fig. 

1.  The  electric  field  of  this  transient  plane  wave  is  assumed  to 
be 

1^(z,t)  *  i  F(t-zcos0/c)  (1) 

where  polarization  is  specified  by  constant  vector  £  and  F(t)  is 
an  unknown  waveform  function  to  be  determined  by  the  requirement  that 
E^(z,t)  excite  a  single-mode,  scattered  field.  The  component  of 
ls1(z,t)  tangential  to  the  wire  surface,  in  the  Laplace  transform 
domain,  is 

^tan(2*s)  "  sine  F(s>e‘SZ  COS8/C  (2) 

This  electric  field  excites  a  transient  induced  current  on  the 
wire,  and  the  induced  current  subsequently  generates  a  transient 
backscattered  electric  field.  The  goal  is  to  synthesize  an  aspect- 
independent  waveform  F(t)  for  the  incident  radar  signal  in  such  a 
way  that  the  backscattered  field,  or  the  return  radar  signal,  from 
the  wire  oriented  at  any  aspect  angle  8  contains  only  a  single 
natural  mode  of  the  target . 

3.  Induced  Current  and  Backscattered  Field 

The  induced  current,  I(z,s),  on  the  wire  can  be  found  as  the 
solution  to  Pocklington ' s  integral  equation  [10] 

Jq  r(z,z' ,s)I(z' ,s)dz'  =  S(z, s)  (3) 

with  the  thin-wire  kernel 


r(z,z' ,s) 


(—  - *2) 


e~sR/c  ,  R  =  V(z-z’)2  +  a2 
4ttR 


(4) 
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and  the  excitation  forcing  function 

S(z , s)  =  -cQ  sind[s  F(s)  e'sz  cos9/c],  (5) 

If  the  only  singularities  of  I(z,s)  in  the  finite,  complex 
s-plane  are  simple  poles  at  natural  frequencies  sa  =  °a  + 
and  it  is  assumed  that  no  entire-function  contribution  exists, 
then  the  SEM  representation  [9,12]  for  the  induced-current  solu¬ 
tion  to  eq.  (3)  is 
N 

I(z,s)  =£  a  (s)  v  (z)  (s-s  (6) 

a=l  a  a  a 

where  v  (z)  is  the  current  distribution  of  the  ath  natural- 
resonance  mode.  Coupling  coefficients  a  (s)  are  obtained  [12] 
from  the  inverse  operator  to  integral  equation  (3)  with  the 
"class-2"  [9]  representation 


A  development  of  this  result  specific  to  the  thin-wire  scatterer 
is  presented  in  the  Appendix;  it  demonstrates  the  apparent  ap¬ 
proximate  nature  of  that  representation  for  aa(s).  Since  the 
above  coupling  coefficient  is  frequency  dependent,  it  leads  as 
necessary  to  an  early-time  response  having  a  forced  component  which 
differs  from  a  pure  sum  of  natural  modes.  If  va(z)  i-s  represented 
by  the  well  known  approximation 


a  (s)  can  be  evaluated  in  closed  form,  and  the  induced  current 
I?z, s)  is  determined  to  be 


-  ** 
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I(z,s)  = 


4it  sin0 


yQL  [2  log(L/a)-l] 


s  F (s)  x 


a=l 


f  i  -  (-i)V 

•SL  COS0/c1[ 

1+1 

.  /  QTf  z\ 

e  1  ni - 1 

Is2  cos2e/ci 

TT7T 
+  a  it  /c  Jl 

-  so.<s-sa>  s:(s's;)J 

sin^  L  / 

The  backscattered  electric  field  can  be  determined  from  the 
far-field  expression 


E  (s)  *  -  SsA(s)sin0 


(10) 


where  A(s)  is  the  vector  potential  in  the  backscattered  direction 
maintained  by  I(z,s),  and  it  can  be  expressed  as 

A  (s)  -  ^  e  '3R“/O/0L  IC^,s)e-sz  dz  (11) 

R  in  eq.  (11)  is  the  distance  between  the  wire  and  a  distant  obser- 
vation  point  in  the  cylinder’s  radiation  zone. 

The  final  expression  for  ES(s)  can  then  be  obtained  as 


-sR  /c 

ES(s)  =  K.,  - — F  (s)  H(s,0) 


where 


K,  = 


1  R 


2 

-L  sm  0 


1  ir2 [2  log (L/a)  -  1] 


,  and 


(12) 


(13) 


F (s)  =  waveform  function  of  the  incident  radar  signal 
Transfer  function  H(s,0)  in  eq.  (12)  is  expressed  as 


N 


H(s,< 


1  -  <-dVsT12[ 

1  i  1  1 

2 

a  L 

1  +  (T/ax)2  s23  1 

sc,(s-sa)  s5(s-s;)J 

where  the  incident -wave front  transit  time  is 
T  =  (L/c)cos8. 


(14) 

(15) 
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It  appears  feasible  at  this  point  to  synthesize  F(s)  required 
to  excite  a  single-mode  scattered  field  by  requiring  Es(s)  to  re¬ 
present  such  a  field  in  eq.  (12)  and  subsequently  solving  for 
F(s)  in  terms  of  known  H(s,9).  This  procedure  was  implemented 
successfully  [14]  for  the  special  case  of  normal  incidence.  With 
obliquely-incident  illumination,  it  becomes  difficult,  due  to  the 
complex  frequency  dependence  of  H(s,6),  to  evatuate  f(t)  by  numeri¬ 
cally  inverse  transforming  F(s);  furthermore  f(t)  is  a  very  ill- 
behaved  signal  in  the  early- time  period.  The  latter  problem  arises 
because  the  above  synthesis  procedure  demands  that  Es(t,0)  con¬ 
sist  of  a  single  natural  mode  at  all  post-incidence  times,  includ¬ 
ing  the  early-time  period  where  the  forced  impulse  response  is 
poorly  behaved  and  difficult  to  approximate.  If  the  frequency- 
domain  synthesis  method  is  abandoned,  then  these  difficulties 
are  overcome  by  adopting  a  time-domain  synthesis  procedure  for 
incident  signals  of  finite  duration  T  .  When  Es(t,0)  is  not  con¬ 
strained  in  the  early- time  period  0  <  t  <  Tg  +  2T,  then  the 
forced  impulse  response  is  not  encountered  and  it  is  found  that  a 
single-mode  scattered  field  can  be  excited  in  the  late-time  period 
t  >  Tg  +  2T  by  well-behaved,  aspect -independent  incident  wave¬ 
forms  of  finite  duration. 

4.  Impulse  Response 

The  impulse  response  h(t,0)  of  the  wire  is  first  determined 
by  inverting  the  transfer  function  as 

h(t,0)  =  l  1  [H (s , 8 ) ]  =  £  ha(t,8).  (16) 

a=l 

h(t,0)  is  a  function  of  aspect  angle  9,  and  h  (t,0)  represents 
that  impulse  response  associated  with  the  a'th  natural  mode.  It 
is  found  that  h^Ct.Q)  can  be  obtained  analytically  in  closed  form 
by  inversion  of  Ha(s,0);  it  is  decomposed  into  an  early-time 
forced  component  which  is  nonzero  only  for  0  <  t  <  2T  and  a  natural 
component,  consisting  of  a  pure  natural  mode,  which  exists  for  all 
0  <  t  <  oo.  These  results  are  summarized  below. 
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During  the  early-time  period  0<t^2T  the  forced  response  is 

"o'''6’  ‘  Tffj7  e°0t  [tVl<t'‘“s,2(tn  C°'  “»t 

+  [uctg1(t)+o<1g2(tn  sin  u0t]  ,  (17) 

while  for  the  late-time  period  of  t>2T  the  natural  response  is 
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h  (t,e)  -  [U(t) G  (t)-2(-l)au(t-T) G  (t-T)+u(t-2T)G  (t-2T) 1 

^  0*  Cl 

'2[Bc.r<V£l<t)-Car<Vf2<t)>  sin  bt 
'2IDc.r(V£l|t)','Ecir,sct)£2<tn  003  bt  <23) 


where 


Ga(t)  =  2e  a  [Aar(sa)  cos  wat-Aai(sa)  sin  «at] 


A  (S  )  —  2  2  s  2  /  B  (s  )  =  s" - ^ - 5 — sr 

u  u  r.*  t  c  A  a.  Ot  a  Z  2  .  ,Z,Z 


b3(sa2-b2) 


a  (sa  +  b  )  ‘ 


2a  (sa  +  b‘) 


Ca^Sa^  ~  _  2  2  5 - »  Da(s  )  = 

2ot  8ct<8a  +b  * 


a  (sq  +b  )  ‘ 


E  (s  )  *  _ ,S. _ 

“  2a2(sa2+b2) 


With  B  (s  )  =  Re{B  (s  )},  etc.  The  first  term  of  eq .  (23)  is  a 
ar'  a  a  a'  n 

natural  mode  and  the  second  and  third  terms  represent  a  forced 
response  which  exists  only  during  the  period  of  0<t<2T  where  T=L 
cos0/c.  Thus,  the  impulse  response,  h(t,e),  consists  of  a  natural 
component  which  is  the  sura  of  natural  modes  and  exists  for  all 
times,  and  a  forced  component  which  exists  only  in  the  early-time 
period  of  0<t <2T. 

Impulse-response  waveforms  for  the  thin-cylinder  structure 
are  indicated  in  Figs.  2  and  3.  Fig.  2  displays  the  impulse 
responses  obtained  by  summing  the  first  10  natural  modes  (first 
layer  of  poles)  in  series  (16)  for  60°  and  30°  aspect  angles.  The 
impulse  response  is  observed  to  depend  strongly  upon  the  aspect 
angle  and  to  consist  of  a  rapidly-oscillating  forced  response 
during  the  early-time  period  0<t<2T  followed  by  the  relatively- 
slowly-varying  natural  response  (a  sum  of  pure  natural  modes)  in 
the  late-time  period  2T<t<®.  Convegence  of  the  impulse-response 
series  was  checked  at  60°  aspect  by  including  999  terms  (pole  lo- 
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cations  estimated  by  an  approximate  recursion  method)  with  results 
as  indicated  in  Fig.  3.  The  result  is  observed  to  be  a  smoothing 
of  the  late-time  response,  which  is  well  approximated  by  the  10- 
term  series,  while  the  rapidly-oscillating,  early- time  response 
persists  with  increased  frequency.  Further  convergence  studies  in¬ 
dicated  that  the  second  and  third  layers  of  poles  provide  essen¬ 
tially  no  contribution  to  the  late-time  natural  response,  while 
their  early-time  nature  is  also  rapidly  oscillatory  and  does  not 
result  in  smoothing  the  total  early- time  impulse  response.  Although 
adequate  convergence  of  the  early-time  impulse  response  could  not 
be  achieved  (this  may  be  related  to  a  Gibb's-type  phenomenon),  it 
is  noted  that  convolution  of  this  response  with  any  relatively- 
smooth  signal,  e.g.,  a  step  or  a  Gaussian  pulse,  will  eliminate 
the  rapid  oscillation  leaving  essentially  the  envelope  of  that 
response.  The  latter  results  are  similar  to  step  and  smoothed- 
impulse  responses  obtained  by  other  investigators. 

The  difficulty  encountered  during  initial  attempts  to  synthe¬ 
size  an  incident  waveform  which  excites  a  single-mode  scattered 
field  at  all  post -incidence  times  is  due  to  the  erratic  behavior 
of  the  approximated  forced  impulse  response  during  its  early-time 
period.  If  the  scatter- field  waveform,  given  by  the  convolution 
E  (t,e)  =  f(t)*h(t,e),  is  required  to  consist  of  a  single  natural 
mode  for  all  t>0,  then  required  incident  signal  f(t)  must  be  such 
as  to  smooth  the  early-time  forced  impulse  response  which  is  in¬ 
tersected  in  the  convolution;  the  result  is  a  poorly-behaved  re¬ 
quired  incident  waveform.  To  overcome  this  difficulty,  a  time- 
domain  synthesis  procedure  for  incident  waveforms  of  finite 
duration  T&  is  developed  which  requires  single-natural-mode  scatter 
field  response  in  the  late-time  period  t>Tg+2T  but  does  not  con¬ 
strain  that  response  in  the  early- time  period  0<t<Tg+2T.  It  is 
emphasized  that  this  method  avoids  the  forced  impulse  response  en¬ 
tirely,  such  that  this  slowly- convergent  component  is  never  required. 

5.  Single-Mode  Excitation 

It  was  demonstrated  in  eq.  (23)  that  the  impulse  response  can 
be  expressed  as 
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h(t,8)  =  C(t,9)  +  an(0)e  n  cos  (tunt  +  4>n  ( 0 ) )  (25) 


where 

£(t,9)  =  forced  response  which  exists  only  during  the 
period  of  0<t<2T, 

N  at 

53  a  (9)  e  cos  (a)  t  +  <£-(8))  =  the  sum  of  natural  modes 
n_-i  n  which  exists  for  all  t, 


a  (0)  =  aspect  dependent  amplitude  of  the  nth  natural 
mode, 

<j>  (0)  =  aspect  dependent  phase  angle  of  the  nth  natural 
n  mode , 

and  an+ju)n  =  sn  *  the  nth  natural  frequency.  h(t,0)  is  shown 
graphically  in  fig.  4a.  _ 

If  the  wire  is  illuminated  with  an  incident  signal,  Ei(t), 
of  duration  as  shown  in  fig.  4b,  the  scattered  field,  Es(t), 
should  appear  like  that  shown  in  fig.  4c,  consisting  of  an  ir¬ 
regular  waveform  for  the  early- time  period  0<t£Te+2T  followed  by 
a  pure  single  natural  mode  for  t>Tg+2T  if  the  waveform  of  Ei(t) 
is  properly  synthesized.  This  phenomenon  can  be  shown  mathemat¬ 
ically  as  follows. 

The  scattered  field,  Es(t,0),  can  be  expressed,  based  on 
the  convolution  theorem,  as 

Es(t,0)  =  fe  Ei(t,)h(t-t' ,0)dt’ 

Jo 


if  N  °rt 

S  (f)  U(t-t\e)  +  2^  an(0)e  n 

•  COB  lw„(t-t *)+♦_ (0 )  J  I  dt' 
n  n  i 


(t-f) 
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Kor  the  late  t  itne  period  t >T  +2T,  the  forced  response  term  does 
not  contribute  to  the  integral  because  £(t-t',0)=O  for  0£t'<T 
if  t>Te+2T.  The  property  £(t,0)  *  o  for  t>2T  has  been  used.  The 
scattered  field  in  the  late-time  period  then  becomes 


„s 


E(  t ,  0 )  = 


f: El<t' 


N  0  ( t— t 1 ) 

£  an(0)e  COSfo) 

n=l 


. .  for  t >Te+2T . 


n(t-t*)+^n(0)]  j  dt1 


(26) 


Equation  (26)  can  be  rewritten  as 


N 


O'  t 


E  (t,0)  *  ]£  an(®)  e  n  fAncos(wnt  +  4>n ( © > )  +Br  sin  (ti>nt  +  ^(8))] 

(27) 


where  the  coefficients  A  and  B  are  defined  as 

n  n 


n 


B 


n 


/Te  i  -ort’ 

E1(t')e  n 


cos  w  t' 
n 

sin  oj  t' 
n 


dt'.  (28) 


It  is  important  to  observe  that  and  Bn  are  independent  of  the 
aspect  angle  0.  It  is  possible  to  choose  an  optimal  E*(t)  in 
such  a  way  that  all  the  coefficients  vanish  except  one,-  by  so 
doing,  Es(t,0)  will  consist  of  a  single  natural  mode. 


6.  Required  Incident  Signals  and  Return  Signals 

It  is  possible  to  choose  an  aspect-independent  E*(t)  to  excite 
a  single-mode  Es(t,e).  Consider  initially  the  class  of  incident 
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waveforms  EX(t)  constructed  as  a  linear  combination  of  natural 
modes  (this  choice  is  generalized  in  Sect.  9)  as 

,  N  at 

EX(t)  =£  e  (bffl  cos  u>mt  +  cm  sin  w  t)  (29) 

m=i 

where  s  *  a  +  ju  is  the  m'th  natural  frequency,  and  b  and  c 
m  m  J  m  n  J  m  m 

are  unknown  coefficients  to  be  determined  based  on  the  requirement 
that  only  a  single-mode  ES(t,6)  be  excited. 

Substituting  representation  (29)  in  eq .  (28)  leads  to 


A 

n 

N 

-L 

mL 

bm 

N  2 

+  E  Mnm  cm 

(30) 

m=l 

m=l 

N 

N 

B 

n 

-E 

Mnm 

bm 

+  E  Mnmcm 

(31) 

m=l 

m*l 

where 


cos 

U) 

t' 

cos 

0) 

n 

m 

cos 

0) 

t’ 

sin 

0) 

n 

m 

sin 

U) 

t ' 

cos 

(1) 

i 

n 

m 

sin 

w 

t' 

sin 

U) 

i 

n 

m 

(32) 


It  is  observed  that  the  ' s  are  explicit  functions  of  incident 

radar  pulse  duration  Ta,  and  T  is  a  parameter  of  freedom  which 

®  i 

can  be  varied  to  obtain  a  desirable  E  (t)  waveform. 

The  unknown  coefficients  b„  and  c  can  be  solved  for  from 

m  m 

eqs.  (30)  and  (31)  as 
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b.  A. 

In  eq.  (33),  matrix  [M^]  is  of  2Nx2N  order,  while  [  c  ]  and  [  g1  ] 
are  two  column  2N  matrices.  i  i 

To  obtain  a  single-mode  scattered  field  (e.g.  the  j'th  mode), 
it  is  required  that  Bj=l  and  Bn*0  for  n#j  and  An==0  for  all  n.  bm 
and  cm  are  easily  determined  from  eq.  (33),  and  E1(t)  is  subsequently 
obtained  from  representation  (29).  With  this  Er(t),  the  scattered 
field,  Es(t,0),  will  be 

s  0it 

E  (t,0)  =  aj(0)e  J  sin(a)jt  +  ( 0 ) )  (34) 

It  is  noted  that  with  this  synthesized  Ex(t),  the  scattered  field, 
s  f 

E  (t,0),  remains  single-mode  for  any  aspect  angle  0,  even  though 

the  amplitude  aj(9)  and  the  phase  angle  <|>j(0)  vary  with  0.  In 

other  words,  when  this  synthesized  E^ft)  illuminates  the  wire,  the 

return  signal  contains  only  a  single  natural  mode  for  any  aspect 

angle  of  the  wire. 

Numerical  results  computed  according  to  this  scheme  are  shown 
in  figs.  5  to  10.  It  was  found  that  when  the  duration  of  the  in¬ 
cident  signal,  Te>  is  one  period  of  the  first  natural  mode,  the 
required  incident  signal  also  becomes  predominantly  single-mode. 

That  is 
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m  1  2  7T  2  IT  »  i  ^  /  L , 

Te  “  T1  =  =  =  oT52Trric7M  =  2,16  (c) 

leads  to  a  desirable  waveform  for  the  incident  signal.  Fig.  5 
shows  the  required  waveforms  for  the  incident  radar  signal  to 
excite  a  return  radar  signal  from  an  arbitrarily  oriented  wire 
containing  only  the  first,  the  second  and  the  third  natural  mode, 
respectively.  For  the  first" mode  excitation,  the  required  incident 
signal  consists  of  one  cycle  of  the  first  natural  mode  and  some 
small  components  of  higher-order  modes  as  shown  in  fig.  5a.  For 
the  second-mode  excitation,  the  required  incident  signal  consists 
of  about  two  cycles  of  the  second  natural  mode  and  some  small  com¬ 
ponents  of  other  modes  as  shown  in  fig.  5b.  This  required  wave¬ 
form  looks  like  a  distorted  second  natural  mode.  For  the  third- 
mode  excitation,  the  required  incident  signal  consists  mainly  of 
about  three  cycles  of  the  third  natural  mode  as  shown  in  fig.  5c. 

It  is  emphasized  that  these  required  waveforms  are  independent  of 

aspect  angle, or  valid  for  any_ aspect  angle  of  the  wire. 

Fig.  6  shows  the  return  radar  signals  from  a  wire  oriented  at 

various  angles,  0=15°,  45°,  60°  and  89°,  when  the  wire  is  il¬ 
luminated  by  the  incident  radar  signal  of  fig.  5a  which  is  synthe¬ 
sized  for  the  first-mode  excitation.  It  is  observed  that  the 
return  signal  for  each  case  of  aspect  angle  remains  that  of  the 
first  natural  mode,  even  though  the  amplitude  and  the  phase  angle 
vary  with  the  aspect  angle.  Also  the  late-time  period  starts  at 
different  time  for  different  aspect  angles.  These  changes  in  the 
amplitude,  the  phase  angle,  and  the  starting  point  of  the  late¬ 
time  period  with  the  change  in  the  aspect  angle  are  not  important 
in  the  practical  detection  of  the  return  signal;  the  most  important 
and  desirable  feature  of  the  return  signal  is  for  it  to  remain 
single-mode.  A  single-mode,  return  signal  can  be  easily  identified 
if  it  is  displayed  visually.  It  is  noted  that  the  return  signal 
during  the  early-time  period  changes  irregularly  with  the  aspect 
angle,  and  it  is  not  shown  in  the  figure  for  the  sake  of  clarity. 

Fig.  7  shows  the  return  radar  signals  from  a  wire  oriented 
at  various  angles  of  15°,  45°,  60°  and  89°,  when  the  wire  is 
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illuminated  by  the  incident  radar  signal  of  Fig.  5c  which  is 
synthesized  for  the  third-mode  excitation.  Again,  the  return 
signal  remains  that  of  the  third  natural  mode  for  any  aspect  angle. 

7.  Target  Discrimination 

When  the  synthesized  incident  signal  for  exciting  a  particular 
natural  mode  of  a  particular  target  is  applied  to  a  different 
target,  the  return  signal  from  the  wrong  target  is  expected  to 
be  significantly  different  from  that  of  a  single  natural  mode. 

Thus,  the  wrongtarget  can  be  sensitively  discriminated.  Two 
numerical  examples  are  given. 

Fig.  9  shows  the  return  radar  signals  from  three  targets,  the 
right  target  (wire) ,  a  wire  5%  longer  than  the  right  target  and 
a  wire  20%  longer,  when  they  are  illuminated  at  30°  aspect  angle 
by  the  incident  radar  signal  of  fig.  5a  which  is  synthesized  for 
exciting  the  first  natural  mode  of  the  right  target.  It  is  observed 
in  fig.  9  that  the  return  radar  signal  from  the  right  target  is 
a  pure  first  natural  mode;  that  from  the  5%  longer  target  displays 
a  slightly  distorted  waveform  and  a  shifted  frequency  from  that 
of  the  first  natural  mode  of  the  right  target;  and  that  from  the 
20%  longer  target  shows  an  irregular  waveform.  Based  on  these 
return  signals,  it  is  easy  to  discriminate  the  wrong  targets  from 
the  right  target. 

Figure  10  shows  the  return  radar  signals  from  the  same  three 
targets  of  fig.  9  when  they  are  illuminated  at  60°  aspect  angle 
by  the  incident  radar  signal  of  fig.  5c  which  is  synthesized  for 
exciting  the  third  natural  mode  of  the  right  target.  The  return 
radar  signal  from  the  right  target  shows  a  pure  third  natural 
mode,  that  from  the  5%  longer  target  displays  an  irregular  ampli¬ 
tude  variation  and  a  shifted  frequency,  and  that  from  the  20% 
longer  target  shows  an  irregular  waveform.  Again,  the  wrong  tar¬ 
gets  can  be  easily  discriminated  from  the  right  target. 


8.  Pulse  Duration  and  Waveform  of  Required  Incident  Signal 
In  Section  6,  it  was  indicated  that  the  waveform  of  the 
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required  incident  signal  for  exciting  a  single-mode  return 
signal  can  be  controlled  by  adjusting  the  pulse  duration, 

T  ,  of  the  incident  signal.  Through  extensive  numerical  cal¬ 
culation,  it  was  found  that  when  Tg  is  equal  to  or  greater 
than  T-p  the  period  of  the  first  natural  mode,  the  aspect- 
independent,  required  incident  sigtial  for  exciting  a  single¬ 
mode  return  signal  consists  mainly  of  that  wanted  natural 
mode,  as  shown  in  fig.  5.  However,  when  Te  is  less  than  , 
the  required  incident  signal  starts  to  have  an  irregular  wave¬ 
form  and  a  high  amplitude.  This  phenomenon  can  be  observed 
in  fig.  11,  where  six  waveforms  for  the  aspect-independent, 
required  incident  signal  for  exciting  the  first  natural-mode 
return  signal  from  an  arbitrarily  oriented  wire  are  shown 
for  signals  possessing  six  different  durations,  Tg  =  0.25  T^, 

0.5  T1#  0.75  ,  T1#  1.5  T.^ ,  and  2  .  It  is  observed  that  for 

the  first  three  cases,  Tq<  the  waveforms  of  the  required 

incident  signal  are  radidly  oscillatory  and  of  high  amplitude, 
while  the  last  three  cases  Tg  >  T^,  have  quite  realizable  wave¬ 
forms  of  mainly  single  mode.  It  is  noted  that  to  eliminate  any 
possibility  of  creating  irregular  waveforms  for  the  cases  of 
Te  <  T7  due  to  numerical  error,  the  required  incident  signal  for 
the  case  of  Tg  -  0.5  T^  was  convolved  with  the  impulse  response 
of  the  wire,  and  the  numerical  result  on  the  return  signal 
turned  out  to  be  a  pure  first  natural  mode  of  the  wire. 

With  this  example,  it  is  easy  to  conclude  that  the  optimum 
pulse  duration  for  the  incident  signal  should  be  a  period  of  the 
first  natural  mode  or  longer.  However,  it  is  found  that  with 
an  incident  signal  of  longer  pulse  duration  the  sensitivity  of 
discrimination  between  different  targets  decreases.  This  is 
easily  visualized  because  as  the  pulse  duration  of  the  incident 
signal  is  increased,  the  situation  approaches  to  the  case  of  the 
continuous  wave  excitation  of  the  target  and  the  return  signal 
will  contain  only  the  excitation  frequency. 

Therefore,  in  the  selection  of  an  optimum  pulse  duration  for 
the  incident  signal  the  following  three  factors  should  be  considered 
(1)  the  pulse  duration  should  be  equal  to  or  longer  than  a  period 
of  the  first  natural  mode  of  the  target,  (2)  an  incident  signal  of 
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longer  pulse  duration  may  be  implemented  with  less  difficulty, 
and  (3)  the  sensitivity  of  target  discrimination  decreases  with 
the  increase  in  the  pulse  duration. 

9.  Uniqueness  of  Synthesized  Waveforms 

The  uniqueness  of  incident  field  waveform  Ei(t),  syn¬ 
thesized  to  excite  a  backscattered  field  consisting  of  a 
single  natural  mode,  is  considered  here.  It  was  demonstrated 
in  the  last  section  that  depends  strongly  upon  its  dura¬ 
tion  Te>  so  the  question  of  uniqueness  must  be  considered 
for  that  class  of  waveforms  with  Tg  specified.  Numerical  re¬ 
sults  in  Sections  6-8  are  all  based  upon  the  natural -mode 
expansion  (29)  for  E1(t)  with  N*10  (N  -  number  of  terms  re¬ 
tained  in  impulse-response  series).  Eq.  (28)  provides  2N 
linear,  algebraic  equations  for  the  2N  amplitude  coefficients 
in  the  series  representation  of  E*.  For  finite  N,  this  solu¬ 
tion  requires  the  backscatter  field  to  consist  of  the  single 
n=j  mode,  while  those  modes  having  l<n<j ,  j<n<N  are  not  ex¬ 
cited.  Modes  having  n>N  are  not  constrained;  this  tranca- 
tion  is  justified  by  the  negligible  contribution  of  such  modes 
to  the  late-time  impulse  response.  The  question  naturally 
arises  whether  a  different  choice  of  basis  functions  in  the 
expansion  for  E*(t)  will  lead  to  the  same  synthesized  wave¬ 
form.  It  is  conjectured  that  E^  is  unique  among  the  class 
of  waveforms  expanded  in  a  complete  basis  set  (the  natural 
modes  are  believed  to  be  complete,  although  they  are  not 
orthogonal)  for  N  infinite.  In  practice,  N  can  be  truncated 
at  some  finite  value  if  the  series  converge  adequately  rap¬ 
idly. 

The  dependence  of  E^  upon  basis  functions  in  its  expan¬ 
sion  with  N"10  was  studied  using  rectangular  pulse,  impulse, 
and  pure  sinusoidal  basis  functions.  Figure  12  indicates  re¬ 
sults  for  E*"(t)  synthesized  to  excite  a  single  mode  (n«l) 
backscatter  field  and  constructed  from  the  various  basis 

functions.  The  signal  duration  was  chosen  as  T  -  T-i ,  and 

i  ex 

it  is  observed  that  E  obtained  with  each  alternative  basis 
set  is  similar  to  the  result  of  the  natural -mode  expansion. 
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Each  of  these  Ex(t)  waveforms  will  excite  the  same  n=l  mode 
backscattered  field,  but  will  not  excite  those  modes  with 
l<n^lO.  These  E1  waveforms  may  excite  different  higher- 
order  modes  with  n>10,  which  have  not  been  constrained.  As 
N  is  increased,  the  E^(t)  waveforms  obtained  with  different 
basis  sets  become  more  nearly  equal,  and  in  the  limit  N-*00 
provide  a  unique  representation  for  E1. 

10.  Conclusion 

It  has  been  demonstrated  that  an  aspect-independent, 
optimal  incident  radar  waveform  E1  of  finite  duration  Te  can 
be  synthesized  to  excite  a  thin-cylinder-target  backscatter 
field  which  consists  of  a  single  natural  mode  of  that  target 
in  the  late-time  period  t>Tg+2T  (T  =  target  transit  time). 

By  constraining  only  the  late-time  target  response,  a  time- 
domain  synthesis  technique  was  developed  which  does  not  re¬ 
quire  knowledge  of  the  forced,  early- time  impulse  response. 

The  optimal  signal  duration  was  found  to  be  near  Tg  =  T^ 

(T-,  =  period  of  first  target  natural -resonance  mode) ,  in 
which  case  the  E  waveform  is  very  nearly  equal  to  that  of 
the  desired  single-natural-mode  return  signal.  Incident 
signals  of  shorter  duration  become  poorly  behaved  with  rapid, 
high-amplilade  oscillations,  while  long  duration  signals 
will  result  in  loss  of  target  resolution  ability.  It  was 
demonstrated  that  a  target-identification  scheme  based  upon 
illuminating  the  target  with  a  waveform  synthesized  to  excite 
a  single  natural -resonance  mode  backscatter  is  capable  of 
sensitive  target  discrimination,  since  the  response  of  a 
wrong  target  with  5%  length  deviation  differs  identifiably 
from  a  single-mode  signal .  Since  the  narrow-band  spectral 
content  of  optimal  incident  waveforms  very  nearly  overlaps 
that  of  the  response  they  excite,  then  this  target  identifi¬ 
cation  technique  should  possess  inherently  good  signal-to- 
noise  ratio  characteristics,  as  well  as  enabling  the  use  of 
narrow-band  filters  and  amplifiers.  Optimal,  synthesized  in¬ 
cident  waveforms  of  specified,  finite  duration  are  unique  when 
represented  by  expansions  in  complete  basis  sets  in  the  limit 
where  all  terms  in  the  impulse -response  series  are  included. 
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The  L.H.S.  of  eq.  (A8)  appears  to  have  sinqulariti.es  in  the  s- 
plane  at  s^,  Sg....,  while  the  R.H.S.  of  the  equation  is  analytic. 

To  clear  this  difficulty,  let's  expand  MOQ (s)  and  (s)  into 

op  p  ot 

Taylor  series'  around  the  respective  poles. 


M88(s) 


>  -  ks(s)l  +  (S-V|w  Mse<s) 

L  J  s-s8  L  J  s= 


+  J(s‘s8,2[^  Mee(s,]s=sB +‘ 


(A9) 


where 


[M80(S3s=s6  / 


vg(z)  dz 


/ 


r(z,z',Sg)  Vg(z’)  dz 


=  0  (A10) 


because  if  Vg(z')  is  the  8th  natural  mode,  no  excitation  is 
needed  to  excite  it,  and  Vg(z')  should  satisfy  eq.  (3)  with 
S (z ,s) =0  . 

[^M88(s,j  s,Sb  -/ 


Vg(z')  dz’i 


Vg(z)  dz 


f  L r~2  r(z,z's)j 

J  o  L3s^  Js=£ 


(All)  _ 
Vg(z')  dz  'j 


Thus,  eq.  (A9)  becomes 


Men15'  *  <S-Ss'  35m88<s)i 


(A12) 


^  • 


(A13) 
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By  expanding  M„  (s)  around  s=s  ,  we  can  similarly  show  that 

ca  ot 


Vls)  *  (s' 


'V  [si  Mea(3)]s=3a+  \  ‘s-sc.>2Uj 


If  we  assume  that 


Vz'>  ’  3in  ' 

we  can  show  after  a  tedious  integration  that 

/IJL  r(z,z',s)  v  (z*)  dz'  =KV  (z) 

o  [as  J  s-sa  a  a 

where  K  is  a  constant.  We  can  then  show  that 

[w  •W‘j  ‘V*»  4/^ 


=  k/L 


vQ (z)  v  (z)  dz  =  0  for  a  *  & 
b  a 


Similarly,  we  can  approximate 


fc?  M«“isisa  4  ° 


for  a  *  8 


Thus , 


V(s)  4  0 


if  a#S 

Equation  (A8)  can  now  be  simplified  to 


as<s,{[a2  MB6  (sjja=Se  2(S-Ss’[^7  Mse(s)]s=s  + 


=  S 


(A14) 

(A15) 

(A16 ) 

dz ' 

(A17) 

(A18) 

(A19) 

,  (s) 


(A20) 


III.IIMIIUM 


a 


The  first  term  of  the  L.H.S.  of  eq.  (A 20)  usually  dominates, 
especially  near  s=Sg.  Therefore,  a^(s)  can  be  approximately 
determined  as 


Vs)  ■ 


S8(s) 


./ 


S (z , s)  V  (z)  dz 
O  P 


[3s  M8B(S)]  s=s6  j^VZ>  dZ  r(z'z’'s]s. 


vR(z': 

SB  6 


(A21) 


which  is  the  result  of  eq.  (7) . 


dz 
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flg.  1.  A  thin  wire  la  Illuminated  by  a  radar 
signal  at  an  oblique  angle. 
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from  a  wire  with  30  and  60  aspect  angles 


The  impulse  responses 

a/a  -  200). 


Waveforms  of  the  Impulse  response  of  a  wire,  the  required 
Incident  signal  and  the  return  signal  with  a  single  natural 
mode  in  the  late-time  period . 


relative  amplitude 


required  waveform  for  the 
first  mode  excitation 


t/ (L/C) 


fig.  5.  Required  waveforms  for  the  incident  radar  signal  to  excite  a  return 
radar  signal  from  an  arbitrarily  oriented  wire  containing  only  the 
the  second  and  the  third  natural  mode,  respectively. 
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fig.  6.  Return  radar  signals  from  a  wire  oriented  at  various  angles, 

8*15°,  45°,  60°  and  89°,  when  it  is  illuminated  by  the  incident 
radar  signal  which  is  synthesized  for  the  first  mode  excitation. 
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4  from  a  wire  oriented  at  various  angles 
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fig.  8. 


Return  radar  signals  from  a  wire  oriented  at  various  angles 
of  I? ,  45° »  60°  and  89° ,  when  it  is  illuminated  by  the  incident 
radar  signal  which  is  synthesized  for  the  third  mode  excitation. 
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ftg.  9.  Return  radar  signals  from  three  targets,  the  right  target  (wire),  a 

wire  5%  longer  than  the  right  target  and  a  wire  20%  longer,  when  they 
are  illuminated  at  30°  aspect  angle  by  the  incident  radar  signal 
which  is  synthesized  for  exciting  the  first  natural  mode  of  the  right 
target . 
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fig.  10.  Return  radar  signals  from  three  targets,  the  right  target  (wire),  a 
wire  5%  longer  than  the  right  target  and  a  wire  20%  longer,  when 
they  are  illuminated  at  60°  aspect  angle  by  the  incident  radar  signal 
which  is  synthesized  for  exciting  the  third  natural  mode  of  the 
right  target. 
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Required  waveforms  for  the  incident  radar  signal  with  various  pulse  durations  to 
excite  an  arbitrarily  oriented  wire  to  produce  a  return  signal  which  contains  only 
the  first  natural  mode  of  the  wire.  These  required  incident  signals  are  aspect- 
independent  . 


Relative  Amplitude 


Natural  modes 


Pure  sinusoids 


Pulse  functions 
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12.  Required  incident  radar  signals  for  exciting  the  first 
natural  mode  backscatter  from  an  arbitrarily  oriented 
wire  of  length  L.  These  signals  are  constructed  with 
natural  modes  of  the  wire,  pure  sinusoids  with  natural 
frequencies  of  the  wire,  pulse  functions  or  impulse 
functions.  The  pulse  duration  of  the  signal  is  equal 
to  a  period  of  the  first  natural  mode  of  the  wire. 


PART  2 


RADAR  WAVEFORM  SYNTHESIS  FOR  EXCITING  SINGLE¬ 
MODE  BACKS CATTERS  FROM  A  SPHERE 


1.  Introduction 

Techniques  of  using  a  short  radar  pulse  to  identify  the  radar 
target  have  been  studied  by  a  number  of  workers  £l-6^.  Typical 
scheme  consists  of  illuminating  the  target  with  a  radar  pulse  and 
then  identify  the  target  by  identifying  it  with  its  natural  modes 
that  are  extracted  from  the  return  signal.  One  of  the  problems 
associated  with  this  scheme  is  difficulty  in  obtaining  accurate 
natural  frequencies  of  the  target  from  a  noisy  return  signal.  In 
a  recent  study  by  our  group,  an  inverse  scheme,  called  the  radar 
waveform  synthesis  method,  has  been  investigated  [7J.  Instead  of 
analyzing  the  return  signal  from  the  target  in  terms  of  its 
natural  mode,  this  new  scheme  synthesizes  the  waveform  of  the  inci¬ 
dent  radar  signal  in  such  a  way  that,  when  it  excites  the  target, 
the  return  radar  signal  contains  only  a  single  natural  mode  of  the 
target.  It  can  be  shown  that  when  the  incident  radar  signal  syn¬ 
thesized  to  excite  a  particular  natural  mode  of  a  preselected 
target  is  applied  to  a  different  target,  the  return  signal  will  be 
significantly  different  from  that  of  the  expected  natural  mode. 

The  wrong  target  can  thus  be  sensitively  discriminated. 

In  this  paper,  the  geometry  of  a  perfectly  conducting  sphere 
is  used  as  the  radar  target.  We  aim  to  synthesize  suitable  inci¬ 
dent  radar  signals  which  can  be  used  to  excite  various  single-mode 
backscatters  from  the  sphere.  To  solve  the  problem,  the  scattered 
field  from  the  sphere  excited  by  an  incident  signal  with  an  arbi¬ 
trary  waveform  is  determined  in  the  spectural  domain.  The  back- 
scattered  field  is  then  determined.  An  approximate  impulse  res¬ 
ponse  in  the  time  domain  is  obtained  by  summing  up  an  infinite 
number  of  natural  modes  based  on  the  Singularity  Expansion  Method. 
That  impulse  response  is  further  approximated  by  the  sum  of  finite 
numbers  of  damped  sinusoids  yielding  an  accurate  result  in  the 
late-time  period. 

It  is  then  demonstrated  that  a  proper  incident  signal  can  be 
synthesized  in  such  a  way  that  when  it  excites  the  sphere,  or  when 
it  convulutes  with  the  impulse  response,  the  return  signal  contains 
only  a  single  natural  mode  of  the  sphere  in  the  late-time  period. 
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A  method  of  synthesizing  required  incident  signals  for  exciting 
various  single-mode  backscatters  is  presented.  Numerical  examples 
are  given  to  show  the  required  incident  signals  for  various  single¬ 
mode  excitations  and  the  resulting  return  signals  which  exhibit 
single  natural  mode  of  the  sphere. 

When  an  incident  signal  synthesized  to  excite  a  particular 
natural  mode  of  a  sphere  is  applied  to  a  wrong  sphere  with  a 
slightly  different  radius  from  that  of  the  right  sphere,  the  re¬ 
turn  signal  from  the  wrong  sphere  is  found  to  be  significantly 
different  from  that  of  a  natural  mode  of  a  sphere.  The  wrong 
sphere  is,  thus,  sensitively  discriminated  from  the  right  sphere. 
This  indicates  the  applicability  of  the  radar  waveform  synthesis 
method  for  radar  target  discrimination. 

2.  Theory 

The  geometry  of  the  problem  is  shown  in  fig.  1  where  a  radar 
signal  propagating  in  the  +z-direction  is  incident  upon  a  per¬ 
fectly  conducting  sphere  which  has  radius  a  and  has  its  center 
located  at  the  origin  of  the  coordinates.  The  electric  field  of 
the  incident  radar  signal  is  assumed  to  be 

Bht,  t)  =  xF(t  -  u  (t  -  (1) 

where  F(t)  is  an  unknown  waveform  function  to  be  synthesized  in 
such  a  way  that  ExCr,t)  excites  a  single-mode  backscatter  from  the 
sphere.  The  Laplace  transform  of  eq.  (1)  can  be  expressed  as 

=  F(s)e"Ya  i  e'YR  c°se  (2) 

A 

where  Y  =  s/c.  When  the  unit  vector  x  is  written  in  terms  of 
spherical  vector  components,  eq.  (2)  becomes 

E*Cr»s)  »  F(s)e’Ya  £sin0  cos<f>  R  +  cos©  cos<{>  0  -  sin<f>|]e  cOS0 

(3) 

Since  E  (r,s)  satisfies  the  wave  equation, 

V  x  V  x  E  +  y2  H  =  0,  (4) 
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^i(r,s)  can  be  expressed  in  terms  of  spherical  vector  wave  func¬ 
tions  (see  Appendix)  as  follows: 


^(r.s)  =  F (s) e” Ya 


(a  +  b 

'•no  In  n  Ain-' 


n  eln; 


(5) 


where 


oln 


sTn&  in(YR)Pn(cose)cos<J,0 


in(YR)- 


■d?1  (cose) 


n 


-JF 


sinij>$ 


rf(i) 

eln 


in(*R)Pn(cos9)cos<}>  R 


SP1 (cos0) 


c  os  4»  6 


+  W  M  tR  inCYR)3 
‘  wlliTe  A  [R  VYR)]  pJi(cos8)sin<j.i 


(6) 


(7) 


The  superscript  (i)  for  Ff  and  if  functions  stands  for  the  use  of 

i  (YR)  which  is  the  first  kind  of  the  modified  spherical  Bessel 
11  1 

function  of  order  n.  Pn  (cos 8)  is  the  associated  Legendre  function 

of  order  n  and  degree  1.  These  functions  are  detailed  in  Appendix. 

The  coefficients  a  and  b  in  eq .  (5)  can  be  determined  by 

11  11  YRcosfi 

comparing  eqs.  (5)  and  (3)  after  e  in  eq.  (3)  is  expanded 

into  spherical  harmonics  ^8^  : 


.n  2n+l 


(-1)  HTn^T  *  br 


(-1) 


n+1 


2n+l 

n(n+l) 


(8) 


Therefore,  the  incident  electric  field  ^("r.s)  can  be  ex¬ 
pressed  as 


£*(r, s)  -  F(s)e 


"  £  t-1'"  K&TT  [«o)n  -  *l\l] 


(9) 


Similarly,  the  scattered  electric  field  from  the  sphere 
^(^r.s)  can  be  expressed  in  terms  of  spherical  vector  functions  as 
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ts  (r,s) 


l:Cs)o 


( -  i ) n  2n+_l__ 
n(n+l) 


where  $2^2  and  n2£2  have  the  same  expressions  of  $^22  and 

om  e±n  r  o^n  eln 

as  given  in  eqs.  (6)  and  (7)  except  with  their  in(YR)  function 

replaced  by  kn(YR)  function  which  is  the  second  kind  of  modified 

spherical  Bessel  function.  kn(YR)  function  is  needed  to  give 

attenuating  behavior  of  the  scattered  field  as  R  approaches  to 

infinity. 

The  coefficient  cn  and  dn  in  eq.  (10)  can  be  determined 
based  on  the  boundary  conditions  on  the  spherical  surface 


Eg  +  Eg  =  0  and  +  E^  =  0  on  R  =  a, 

to  be 

c  -  >^a) 

n  kn(Ya) 


a 

W 

R  in(YR)= 

R=a 

a 

W 

R  kn(YRX 

R=a 

(ID 

(12) 


l  , 


The  final  expression  for  the  scattered  electric 


field  is 


r,s) 


F(s)e 


-  Ya 


(-D 


n  2n+l 
n (n+1 ) 
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'00 
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^j(R  kn(VR)) 


(13) 


The  backscattered  electric  field  in  the  far  zone  of  the 
sphere  should  be  in  the  x-direction,  parallel  to  the  direction  of 
the  incident  electric  field,  and  can  be  found  from 
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E3  -  *  [E3] 


R=oq  ,  e=1T,  <j>  =  1T 


Since  as  R  approaches  to  infinity, 

-\R  - VR 

kn(YR)“"7  W  e  and  W  7k(R  kn(YR})~*  '  7  W  e 


|pl(cos6) | 

L  sin0  J 


=  -  JLlI 


f-nn  r3pJ(cose)l 


=  C-i) 


-n(n+l) , 


as  R -►<*>,  0  =  ir  and  4)=ir  ,  F^k^  and  >f^k^  functions  become 

ass  -  *.S-i  ?  4  e-VR“^(n.1). 

Thus,  the  backscattered  electric  field  becomes 


F(s) 


,e'Y(R“+a)  ,n 

:  *F (s) — 2-i  (2n+l)  - g-p - =\ - 

-  s*  l  va)*Ki'4. 


(2n+l) 


oo  n=j 


V*)  sWlh 


Using  the  following  Wronskin, 


Vx)lr[x  Vx>]  •  [x  Vx)]  =  •  7© 


eq.  (15)  can  be  simplified  to  be 

-  2  *Y(R0o+a)  - 
E5 (s)  =  F(s)  g— - -V  2 


Y  ^  ^[#R  VYR))JR=a 


The  backscattered  electric  field  £s  given  in  eq.  (17)  can  be 
developed  further  because  the  modified  spherical  Bessel  functions 
can  be  expressed  with  truncated  power  serieses  as  follows  [9]. 
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k„<»  -  if 

^  (n+8) ! 

ehn-ej!  (2s) 

(18) 

■h  - 

_  *  (n+  6) !  nJunrS 

7  e  iL  ITTn %T  (1V(2c} 

(19) 

If  eqs.  (18)  and  (19)  are  substituted  in  eq.  (17),  we  have  the 
final  expression  for  E  (s)  as 


^(s) 


''la 

X  I  Rm  e 


Y(R00-a)p(s) 


'f'  (2n+l)c2n 


(20) 


where 


£-K)  ■  k 

(n+B ) !  1  ,n-B 

(21) 

8  t (n-e)  !  2b  4 

g"u)  ■  £ 

Cn+Bli  1  (B+r)rn'e 

,  b!  (n- ejr  2b  v-b 

(22) 

?  =  Ya 

and  Y  =  s/c. 

=  distance  between  the  observation  point  and  the 
center  of  the  sphere. 


The  backscattered  electric  field  can  also  be  expressed  as 


£s(s)  =  -x  -  e"s{R»*a)/c  F(s)H(s) 


where 


H(s)  ’  £  wlth  ‘-Is 


(23) 


(24) 


In  eq.  (23),  F(s)  is  the  unknown  function  describing  the  waveform 
of  the  incident  radar  signal,  and  H(s)  is  the  transfer  function  of 
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the  sphere.  The  transfer  function  H(s)  as  expressed  in  eq.  (24) 
is  a  poorly  converging  infinite  series  and  its  evaluation  for  large 
value  of  c  requires  sophisticated  Watson's  transformation  £loJ  or 
extensive  Fourier-numerical  method  £llj.  To  synthesize  an  incident 
radar  signal  which  excites  a  single-mode  backscatter  from  a  sphere, 
it  is  necessary  to  obtain  a  reasonably  accurate  impulse  response 
which  is  the  Laplace  inverse  transform  of  H(s).  This  is  done  in 
the  next  section. 

3.  Impulse  Response 

The  impulse  response  h(t)  is  obtained  by  inverting  H(s).  In 
the  process  of  inverting  H(s),  the  roots  of  f  (5)  and  gn(c)  func¬ 
tions  in  eq.  (24)  are  computed  based  on  Muller's  algorithm.  There 
are  n  roots  for  each  fn(C)  function  and  n+1  roots  for  each  gn(c) 
function,  as  can  be  seen  easily  from  eqs.  (21)  and  (22).  If  these 
roots  are  plotted  in  the  s-plane  (or  the  £-plane),  they  can  be 
grouped  into  branches  of  roots  as  shown  in  fig.  2.  Except  those 
roots  lying  on  the  negative  real  axis,  all  other  roots  are  in  con¬ 
jugate  pairs.  These  roots  indicate  the  locations  of  the  simple 
poles  of  H(s)  in  the  s-plane.  It  is  noted  that  the  roots  in  one 
branch  come  from  either  fn(0  or  gnCO  but  they  do  not  belong  to 
the  same  index  n.  For  example,  the  first  branch  (and  other  odd- 
numbered  branches)  of  roots  come  from  gn(t)  functions  of  various  n, 
and  the  second  branch  (and  other  even-numbered  branches)  of  roots 
come  from  fn(?)  function  of  various  n.  Through  numerical  calcula¬ 
tion,  we  found  that  this  kind  of  regrouping  the  roots  in  the  s-plane 
provides  an  interesting  information;  asymptotically,  the  roots  be¬ 
longing  to  the  same  branch  have  a  simple  arithmetic  relation  be¬ 
tween  their  locations  on  the  s-plane,  and  the  residues  of  H(s)  at 
these  roots  (or  poles)  possess  a  simple  geometric  relation  between 
their  complex  amplitudes. 

For  simplicity,  we  can  designate  the  roots  with  the  symbol , 

C  ,  where  the  subscript  i  represents  the  ith  branch  and  the  sub¬ 
script  j  means  the  jth  root  of  the  branch.  Also  we  can  use  j=0 
to  designate  those  roots  lying  on  the  negative  real  axis.  In 
general,  we  can  express 
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^ij  ~  aij  +  ^ij  ~  °ij  '  anc*  ^iO  ”  aiO’ 


The  transfer  function  H(s)  contains  only  simple  poles  located 
at  these  roots  and  the  residues  of  H(s)  evaluated  at  these  poles 
are  given  by 


[Res  H(s)  ]at  [Res  H(s)]at 


Some  lower-order  roots  and  the  corresponding  residues  of  H(s)  are 
tabulated  in  Table  1. 

The  transfer  function  H(s)  can  be  expanded  into  an  infinite 
series  as 


H  (s) 


(25) 


The  summation  of  the  second  term  of  eq.  (25)  starts  from  i=2  be¬ 
cause  there  is  no  root  of  the  first  branch  lying  on  the  negative 
real  axis.  The  upper  limit  of  the  summation  over  i,  N,  is  dictated 
by  the  maximum  number  of  the  branches  needed  to  be  considered. 
Numerically  it  was  found  that  a  reasonably  accurate  solution  can  be 
obtained  with  N  less  than  3. 

The  impulse  response  can  then  be  obtained  as 


h(t) 


■ «  £  £ 


2  Re(ai;je  * 


(=)  £  n< 


5iOT 


(26) 


where  x  *  t/(a/c)  is  a  normalized  time,  and  each  term  of  eq.  (26) 
represents  a  natural  mode  of  the  sphere. 

Examining  Table  1,  one  can  observe  a  simple  arithmetic  relation 
between  the  values  of  roots  (or  the  locations  of  poles)  belonging 
to  the  same  branch  if  j  is  bigger  than  10.  That  is, 

=  ASj  ■*  constant,  if  j  >  10.  For  example, 

-  -0.035  +  j0.98  and  Ai^  =  -0.078  +  j0.96  for  j  >  19. 
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Similarly,  there  is  a  simple  geometric  relation  between  the 
residues  at  these  poles  of  the  same  branch  if  j  is  bigger  than  10. 
That  is,  =  R^  constant,  if  j  >  10.  For  example, 

Rx  =  1.09  e-^ 1  ‘ 18  and  R2  =  1.19  e^1*22  for  j  >  19. 

In  view  of  these  simple  relations,  the  infinite  sum  of  the 
modes  in  eq.  (26)  can  be  divided  into  two  parts;  the  first  part  is 
the  sum  of  the  first  10  or  20  terms  and  the  second  part  is  the  sum 
of  the  rest  of  the  terms.  For  example,  we  can  sum  the  modes  coming 
from  the  poles  of  the  first  branch,  i=l,  as 

[I'M]  i  -  (|)2  Re  U  axje  ^  ,  (c)  n,  ^  ,  [;e  ?UT 

(27) 


Now  the  simple  relations  between  roots  and  residues  for  large  j  can 
be  used  to  approximate  the  last  term  of  eq.  (27)  as  an  infinite  geo¬ 
metric  series  which  is  then  summed  up  to  be 


The  fact  that  A^  contains  a  negative  real  part  has  been  used,  and  eq.  (28)  is 
valid  for  r  >  0.  Fort=0,  the  impulse  response  can  be  shown  to  be  infinite, 
directly  from  eq.  (24)  based  on  the  initial  value  theorem. 

If  numerical  values  of  a^  ^  ^g,  and  A^  are  substituted  in  eq.  (28), 
we  have 


2  Re  [£„  a'ie  'ljT]  = 


89.34e 


0 

-1. 797t 


_ Ceos (19. 3 2t  -0.485)-1.09e'0,035Tcos(18.34T -1.66)3 

1  +  1.19e'O,°7T-2.18e'0d38Tcos(0.98T  +  1.18) 


(29) 
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The  first  term  of  eq.  (27)  can  be  easily  summed  up  numerically, 
and  when  it  is  combined  with  eq.  (29)  the  part  of  the  impulse  res¬ 
ponse,  h(t),  which  is  contributed  by  the  modes  of  the  first  branch 
poles,  can  be  computed  as  a  function  of  time.  It  was  found  that  a 
reasonably  good  result  can  be  obtained  just  considering  the  contri¬ 
bution  from  the  poles  of  the  first  branch. 

A  numerical  example  is  given  in  Figs.  3  and  4.  Fig.  3  shows 
the  impulse  response  computed  from  the  sum  of  the  first  19  pairs 
of  the  poles  of  the  first  branch,  or  the  first  term  of  eq.  (27). 

It  shows  a  strong  oscillatory  response  in  the  early-time  period, 
but  no  creeping  wave  peak  is  observed.  This  result  is  wrong, 
judging  from  the  existing  results.  However,  if  the  term  given  in 
eq.  (29),  which  represents  the  contribution  due  to  the  rest  of  the 
poles  of  the  first  branch,  is  added  to  the  contribution  from  the 
first  19  pairs  of  the  poles,  a  surprising  result  is  obtained;  the 
strong  oscillatory  response  during  the  early-time  pariod  is  can¬ 
celled  and  a  sharp  peak  representing  the  creeping  wave  contribution 
appears  at  x  =  t(c/a)  *  5.25,  as  shown  in  Fig.  4.  It  is  noted  that 
an  impulse  at  t=0  is  added  in  Fig.  4,  as  it  should  be.  The  impulse 
response  shown  in  Fig.  4  agrees  with  the  existing  results.  An  ap¬ 
proximate  impulse  response  based  on  Fourier-numerical  method  £llj 
is  included  in  Fig.  4  for  comparison.  If  the  contribution  from 
the  poles  of  other  branches  is  considered,  the  accuracy  of  the  im¬ 
pulse  response  during  the  early-time  period  can  be  improved.  It 
is  noted  that  if  more  than  19  pairs  of  poles  in  the  first  branch 
are  considered  without  the  compensation  of  a  corresponding  correc¬ 
tion  term,  the  oscillatory  response  in  the  early-time  period  be¬ 
comes  stronger  instead  of  weaker.  This  is  due  to  the  fact  that 
the  amplitude  (residue)  of  the  mode  increases  with  the  order  of  the 
mode  (Table  1).  For  many  practical  applications  the  result  of 
Fig.  4  is  sufficient. 

For  our  purpose  of  synthesizing  an  incident  radar  signal  which 
excites  a  single-mode  backscatter  from  a  sphere,  it  is  desirable 
to  obtain  an  approximate  impulse  response  in  the  form  of  a  truncated 
sum  of  natural  modes.  To  do  so,  it  is  necessary  to  approximate  the 


correction  term  of  eq.  (29)  with  two  damped  sinusoids  which  pos¬ 
sess  the  forms  of  natural  modes.  This  step  was  accompl ished  nu¬ 
merically  by  the  cut  and  try  approach.  We  found  the  two  damped 

sinusoids,  1.34  e'2'°T  cos  (18 . 8 9 1 t  +  .4838)  and 
-  55t 

.127  e  ‘  cos  (19.50x  -  .7159),  can  approximate  the  correction 
term  of  eq.  (29)  quite  well  for  x  >  2.  Since  we  aim  to  synthesize 
a  required  incident  radar  signal  for  exciting  a  single-mode  back- 
scatter  in  the  late-time  period,  we  do  not  need  the  information 
on  the  early-time  behavior  of  the  impulse  response.  The  approxi¬ 
mate  impulse  response  constructed  with  the  first  19  natural  modes 
of  the  first  branch  and  two  damped  sinusoids  which  approximate  the 
correction  term  of  eq.  (29)  is  shown  in  Fig.  5.  This  approximate 
impulse  response  approximates  the  true  impulse  response  quite 
accurately  for  x  >  2.  Thus,  it  will  be  used  in  the  synthesis  of 
the  incident  radar  signal  for  exciting  a  single-mode  backscatter 
in  the  late-time  period. 


4.  Excitation  of  Simple-Mode  Backscatter 

The  approximate  impulse  response  of  a  sphere  can  be  represented 
by  the  first  19  natural  modes  from  the  first  branch  of  poles,  two 
damped  sinusoids  which  approximate  the  correction  term  of  eq.  (29) 
for  x  _>  2,  and  a  term  which  compensates  the  error  for  the  early- 

time  period  of  0  <  x  <  2.  Symbolically,  we  can  express  the  impulse 
response  h(x)  as 


h(r) 


■2 


o  x 
n 


n 


cos 


(d)  T 

n 


<p  ) 


e(T) 


(30) 


In  eq.  (30),  the  first  19  terms  of  the  summation, 


2 


an  e 


0nT 


cos  (wnx  +  4>n), 


modes  from  the  first  branch 
summation,  / 

mo 

sinusoids  which  approximate 
x  >_  2 ,  i .  e.  , 


a  x 

a  e  n  cos 
n 


represent  the  sum  of  the  first  19  natural 
of  poles.  The  last  two  terms  of  the 
(«nt  +  4>n)  ,  represent  the  two  damped 
the  correction  term  of  eq .  (29)  for 
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a2Q  e  20  cos((i>2qt  +  4>20)  =  (|)  1.34  o' 2  * 0  ’cos  U  8 . 891  i  +  .48381 

Cf  ^  i  T  q  q 

a21  e  cos(a>21T  +  <j>21)  =  (|)  0.127  e'  *  Tcos  (19. 50t  -  .7159). 


The  term  C(x)  exists  only  during  the  period  of  0  £  x  £  2 ,  and  is 
the  term  to  be  added  to  the  summation  term  of  eq.  (31)  to  yield  an 
accurate  impulse  response  for  that  period  of  time  because  with  only  the 
summation  term  of  eq.  (31)  it  does  not  give  accurate  result  for  h(t) 
during  that  early-time  period.  The  function  £(x)  is  difficult  to 
determine  but  it  is  not  needed  if  we  only  aim  to  produce  a  single¬ 
mode  backscatter  in  the  late-time  period. 

We  now  aim  to  synthesize  an  incident  electric  field  E1(t)  of 
duration  xe  in  such  a  way  that  when  it  illuminates  the  sphere,  the 
backscattered  electric  field  Es(t)  consists  only  of  a  single 
natural  mode  in  the  late-time  period  of  x  >  r  +  2.  The  backseat- 
tered  electric  field  E  (t)  can  be  expressed,  based  on  the  convolu¬ 
tion  theorem,  as 

ES(t)  =  f  6  E*(x')h(x  -  x")dx' 

*^o 


an(x-x") 

e  cos  (u>n(x-x'  ) +4>n)  +  £(x  -  x  ')  dx' 


(31) 


If  the  normalized  observation  time  x  >  x„  +  2,  the  term  e(x-x') 
does  not  contribute  to  the  integral  because  £(x)  =  0  for  x  >  2. 
Thus,  eq.  (31)  becomes 

es(t)  .  jtW(t') 


[£ 


(x  -x'  ) 


an  e 


cos(wn(x-x  ) 


■ 

+  4>  ) 
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where  the  coefficients  A  and  B  are  defined  as 

n  n 


;•}  •  f:  ■■ 


-  0  T  COS  W  T 

( t  )  e  dr 

sin  w  t" 
n 


Based  on  eq.  (33),  it  is  now  possible  to  choose  an  optimal  E1(r) 
in  such  a  way  that  all  the  coefficients  vanish  except  one;  by  so 
doing,  Es(t)  will  consist  of  a  single  natural  mode. 


5.  Required  Incident.  Signals  and  Return  Signals 

It  is  possible  to  choose  an  E1(t)  to  excite  a  single-mode  Es(t) 
Consider  an  incident  electric  field  E^x)  constructed  from  a  linear 
combination  of  21  damped  sinusoids,  19  natural  modes  and  two  damped 
sinusoids  approximating  the  correction  term  of  eq.  (29),  as  that 
appeared  in  eq.  (30): 

21  or 

E1(t)  =  2  e  m  (b  cos  w  x  +  cm  sin  ^r)  (34) 

m=l 

where  c  =  cr  +  jw  is  the  m'th  natural  frequency  and  and 

in  m  in  l  &  x 

are  the  equivalent  values  of  the  two  damped  sinusoids,  and  bm 

and  c_  are  unknown  coefficients  to  be  determined  based  on  the  re- 

m  s 

quirement  that  only  a  single-mode  E  (x)  be  excited. 

Substituting  representation  (34)  in  eq .  (33)  leads  to 


21  21 

=  12  M*  K  +  2  Mn 
*—J  nm  m  n 


nm  cm 


v  -  £  k 


Mnm  bm  +  Mnm  cm 
m= 1  m=l 


where 
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(37) 


It  is  observed  that  the  M1  's  are  explicit  functions  of  incident 

nm  r 

radar  pulse  duration  x  ,  and  x  is  a  parameter  of  freedom  which 

e  e  i 

can  be  varied  to  obtain  a  desirable  E  (x)  waveform. 

The  unknown  coefficients  b  and  c  can  be  solved  for  from 

m  m 

eqs.  (35)  and  (36)  as 


To  obtain  a  single-mode  scattered  field  (e.g.,  the  j'th  mode), 

it  is  required  that  B.  =  l  and  B  =0  for  n?<j  and  A  =0  for  all  n.  b 

1  n  n  i  m 

and  cm  are  easily  determined  from  eq.  (38),  and  E  (x)  is  subse¬ 
quently  obtained  from  representation  (34).  With  this  E1(x),  the 
scattered  field,  Es (x ),  becomes  single-mode  and  can  be  expressed  as 


15- 


a-  t 


ES(t)  =  e  J  sin  (ukt  +  <Jk  ) 


(39) 


Some  numerical  results  on  the  required  incident  signals  for 
exciting  various  single-mode  backscatters  and  resulting  single¬ 
mode  return  signals  are  shown  in  Figs.  6  and  7.  Figure  6  shows 
the  required  incident  signal  for  exciting  the  first-mode  backscatter 
(^1  =  -0.5  +  j0.866)  and  the  return  signal  which  indeed  shows  the 
first  natural  mode  in  the  late-time  period  of  t  >  9.  The  duration 
of  the  required  incident  signal  is  set  to  be  one  period  of  the  first 
natural  mode. 


T 

e 


T  =  1  =  -  2  TT 

1  0 . 866 


7.26. 


The  waveform  of  the  required  incident  signal  is  found  to  contain 
a  rapidly  oscillatory  component  in  the  initial  stage.  However, 
the  return  signal  contains  only  the  much  slower  varying,  first 
natural  mode  in  the  late-time  period.  This  phenomenon  is  different 
from  the  case  of  a  thin  wire  where  the  required  incident  signal 
for  single-mode  excitation  consists  mainly  of  the  wanted  natural 
mode  ([7].  The  reason  is  that  the  natural  modes  of  a  thin  wire  are 
nearly  orthogonal  while  that  of  a  sphere  are  not  orthogonal  due  to 
their  large  damping  coefficients.  It  is  noted  that  the  return 
signal  was  obtained  by  convoluting  the  required  incident  signal 
with  the  approximate  impulse  response  given  in  eq,  (30).  The  early- 
time  part  of  the  return  signal  exhibits  an  irregular  waveform  and 
is  not  shown  in  the  figure  for  the  sake  of  clarity. 

Figure  7  shows  the  required  incident  signal  for  exciting  the 
third-mode  backscatter  (s3  =  -0.843  +  j  2 . 7  58)  and  the  return  signal 
which  contains  only  the  third  natural  mode  in  the  late-time  period 
of  t  >  9.  The  required  incident  signal  has  a  duration  of  one  period 
of  the  first  natural  mode,  and  its  waveform  consists  of  a  rapidly 
oscillatory  component  superimposed  on  a  slowly  varying  component. 

The  third-mode  return  signal  was  created  by  the  numerical  convolu¬ 
tion  of  the  required  incident  signal  with  the  approximate  impulse 
response . 

The  required  incident  signals  for  exciting  the  first  and  the 
third-mode  backscatter  appear  to  be  somewhat  irregular  in  waveform. 
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However ,  they  can  be  constructed  with  21  natural  modes  of  appro¬ 
priate  amplitudes  and  phase  angles  as  shown  in  Table  2.  It 
appears  that  higher  amplitudes  for  higher-order  natural  modes  are 
needed  to  construct  the  required  incident  signal. 

The  relation  between  the  signal  duration  and  the  waveform  of 
the  required  incident  signal  was  also  studied.  It  was  found  that 
when  the  signal  duration  is  shortened,  the  waveform  of  the  required 
incident  signal  resembles  that  shown  in  Fig.  6  or  Fig.  7  but  it 
contains  a  more  rapidly  oscillatory  component  in  the  initial  stage. 

5.  Target  Discrimination 

To  show  the  capability  of  target  discrimination  of  this  method, 
two  numerical  examples  are  given  in  Figs.  8  and  9.  Figure  8  shows 
the  return  signals  from  two  spheres,  the  right  sphere  and  a  wrong 
sphere  which  radius  is  10%  smaller  than  that  of  the  right  sphere, 
when  they  are  illuminated  by  the  required  incident  signal  for  ex¬ 
citing  the  first-mode  backscatter  from  the  right  sphere  as  that 
shown  in  Fig.  6.  It  is  observed  that  the  return  signal  from  the 
right  sphere  is  a  pure  first  natural  mode  while  that  from  the 
wrong  sphere  shows  distortions  in  its  waveform.  Figure  9  shows  the 
return  signals  from  the  same  two  spheres  when  they  are  illuminated 
by  the  required  incident  signal  for  exciting  the  third-mode  back¬ 
scatter  from  the  right  sphere  as  that  shown  in  Fig.  7.  The  return 
signal  from  the  right  sphere  is  a  pure  third  natural  mode  while 
that  from  the  wrong  sphere  exhibits  an  irregular  waveform.  From 
these  two  examples,  it  is  evident  that  the  wrong  target  can  be  sen¬ 
sitively  discriminated  from  the  right  target  if  the  incident  sig¬ 
nals  are  properly  synthesized  for  single-mode  excitation.  It  is 
noted  that  if  the  wrong  target  were  a  thin  wire  or  some  other  non- 
spherical  object,  the  return  signals  from  the  wrong  target  would  be 
entirely  different  from  the  natural  modes  of  a  sphere,  and  the  dis¬ 
crimination  of  the  wrong  target  from  the  right  target  would  be  tri¬ 
vial  in  that  case. 

6.  Conclusion 

It  has  been  demonstrated  that  by  expressing  the  backscattered 
field  of  a  sphere  in  terms  of  its  natural  modes,  it  is  possible  to 
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synthesize  an  incident  radar  signal  of  appropriate  waveform  in 
such  a  way  that  when  it  illuminates  the  sphere,  the  return  signal 
contains  only  a  single  natural  mode  of  the  sphere  in  the  late- time 
period.  When  that  synthesized  incident  signal  is  applied  to  a 
wrong  target,  the  waveform  of  the  return  signal  will  be  signifi¬ 
cantly  different  from  that  of  natural  modes  of  the  sphere;  thus, 
the  wrong  target  can  be  sensitively  discriminated  from  the  right 
target.  Since  a  single-mode  radar  return  contains  a  narrow  fre¬ 
quency  band  it  may  lead  to  improved  signal- to-noise  ratio  and  can 
be  processed  by  narrow-band  filters  and  amplifiers. 
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Appendix:  Solutions  of  Wave  Equation  in  Terms  of  Spherical  Vector 

Wave  Functions. 


Consider  the  wave  equation  for  the  Laplace  transformed  F  field 
VxVxE+Y2E=0  (Al) 


where  Y  =  s/c. 

To  solve  eq.  (Al) ,  we  consider  first  the  corresponsing  scalar  wave 
equation, 

V2g  -  Y2g  =  0  (A2) 


Assuming  g(R,0,<f>)  =  g^  (R)  g2  (9 )  gj  (4>)  »  and  by  the  technique  of  vari¬ 
able  separation,  g ^  CRD  , g 2  ( 0 )  anc*  83  (<fr )  can  be  shown  to  satisfy 
the  following  equations. 


R 


,  2 

d  g] 


+  2R  ^  +  [-Y2R2-n(n+l)]  g2  *  0 


(A3) 


1  d 
sm0  3F 


dg2 


(sin0  gg^)  +  [n(n+l)  - 


2 

m 

~  "~T ' 
s  m  0 


(A4) 


d2&3  2 

+  m  83  =  ° 


(A5) 


The  solution  for  g1(R)  is  given  by 


n+^R> 

Kn+h(YR) 


in(VR) 


kn(YR) 


r\ 


(A6) 


where  In+jj(yR)  is  the  first  kind  of  modified  Bessel  function  of 
order  n+*j,  Kn+Jj(YR)  is  the  second  kind  of  modified  Bessel  function 
of  order  n+*5.  in(YR)  and  kn(YR)  can  be  considered  as  the  first  and 
second  kind  of  modified  spherical  Bessel  functions  of  order  n. 
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The  solution  for  g 2 ( e )  is 

g2(6)  *  P„(cos8)  U”) 

where  Pm(cos0)  is  the  associated  Legendre's  function  of  order  n  and 
degree  m. 

The  solution  for  g3(e)  is 


g*(9)  = 


COS  Tliip' 

=  -  ► 
sin  m<j> 

. 


Therefore,  the  general  solution  for  g(R,0,<j>)  is 


g(R»9  >$)  -  EE 


n  m 


VYR>1  pm 

k  (YR) 
n 


(cose) 


cos  m<}> 
sin  m<j> 


It  is  possible  to  construct  from  g(R,0,<|>)  two  possible  solu¬ 
tions  for  I?  which  satisfy  eq.  (Al)  and  Maxwell’s  equations  as  fol¬ 
lows.  The  first  solution  for  £  is 


Ft  (R,  8 ,  d>)  *  V  x  (£g) 
and  the  second  solution  for  £  is 
(R,9,<t>)  -  y  V  X  ft  (R,0,$) 


(A10) 


(All) 


These  two  solutions  can  be  shown  to  satisfy  eq.  (Al)  by  direct  sub¬ 
stitution. 

The  substitution  of  eq.  (A9)  in  eq .  (A10)  gives 


-  I  5 
omn  " 


i  (YR)1  n. 

n  l  pm 
kn(YR)J  n 


(cose ) 


rin(vR)  1  _pJ 

\VYR>  j 


sin  I 

r  9 

cos  m<j>  I 


(cose)  cos 
515  * 

sin  miji 


(A12) 
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The  substitution  of  eq.  (A12)  in  eq.  (All)  leads  to 


tf  =  nMI  f  pm 

§mn  ™  \kn(YR)  J  n 


(COS0) 


cos  m<J> 


sin  m<{> 


.  i  a  /  R  V™'  1  e>0^ 
(R  kn(rR)  J  96 


TRsinS  3R 


R  in(YR) 
R  kn(YR) 


P™  (COS0) 


cos  m(j> 
sin  m<j> 


sin  m<j> 
cos  m4> 


(A13) 


The  subscript  e  and  o  stand  for  the  even  and  odd  $  and  $  functions. 
The  general  solution  for  £  is  any  combination  of  even  and  odd  Ft 
and  it  functions  of  any  m  and  n. 

It  is  noted  that  these  solutions  deviate  from  the  conventional 
Mie  series  ^8j  in  that  the  present  solutions  use  modified  spherical 
Bessel  functions  while  Mie  series  uses  ordinary  spherical  Bessel 
functions . 
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fig.  2  Natural  frequencies  of  a  conducting  sphere;  the  poles  are 
grouped  into  branches  in  the  S-plane 

Table  1:  Poles  of  the  first  and  second  branches  shown  in  the  above  figure 

and  the  corresponding  residues. 


poles  of  1st 
branch 


residues  at  c 


poles  of  2nd 
branch 
C2j  *  °2j  + 


residues  at  c 


-.9540 
-1.048 
1.129 
-1.201 
-1.267 
-1.327 
-1.382 
-1.434 
-1.483 
-1 .528 
-1.572 
~1 .613 
-1.653 
-1.69 1 
-1.727 
-1.762 


.8660 

I. 807 
2.758 
3.715 
4.676 
5.642 
6.610 
7.580 
8.553 
9.527 
10.50 

II. 48 
12.46 
13.44 
14.42 
15.40 
16.38 
17,36 
18.34 


0. 

-2.757 
-1.879 
3.734 
6.145 
-.8226 
-9.734 
-6.852 
7.700 
15.53 
2.777 
-17.67 
-18.13 
7.004 
28.20 
14.88 
-21 .84 
-35.97 
-3.843 


1.732 

.5773 


WW yA 


-3.825 

2.920 

8.327 

2.529 

-9.726 

-11.51 

3.468 

17,93 

10.41 

-14,09 

-24.65 

-3.078 

27.57 

26.61 

-11.00 

-40.83 


-1.000 
-1.500 
-1 .839 
-2.104 
-2,325 
-2.516 
—2 , 686 
-2.839 
-2.979 
-3.109 
-3.230 
-3.343 
-3.450 
-3.551 
-3.647 
-3.739 
-3.827 
-3.912 
-3,993 
-4,071 


.8660 

I .  754 
2.657 
3.571 
4,493 
5.421 
6.354 
7,292 
8,233 
9.177 
10.12 

II. 07 
12  03 
12.98 
13.94 
14.89 
15,85 
16,81 
17.77 


2  j 


3.000 

0. 

-18.05 

-10.29 

45.70 

58.40 

-57.95 

-163.4 

-15.41 

284 . 3 

250.3 
-278.2 
-627, 
-60.65 

919.4 
828.8 

-697,2 
-1771 . 
-425.9 
2165. 


0. 

8.660 
2.364 
-31 .07 
-27,74 
57.23 
104.2 
--37.51 
-228.3 
-110.1 
309 . 9 

429.4 
-161 .8 
-807.6 
-396.4 

902.5 
1313. 

-259.0 

-2096, 

-.1325, 
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fig.  3  The  impulse  response  of  a 
computed  from  the  first  19 
branch . 


Scattered  Field 


o 


fig.  4  The  impulse  response  of  a  conducting  sphere  of  radius  a  computed 
from  the  first  19  poles  of  the  first  branch  and  the  correction 
term  of  eq.(29).  The  impulse  response  shown  by  the  dotted  line 
is  obtained  with  the  Fourier-numerical  method. 
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fig.  5  Approximate  impulse  response  for  t>2  constructed  with 
the  first  19  natural  inodes  of  the  first  branch  poles 
and  two  damped  sinusoids  which  approximate  the  correction 
term. 


Relative  Amplitude  Relative  Amplitude 
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Table  2.  Natural  Modes  Used  to  Construct  the  Required  Incident 

Signals  for  Exciting  the  First  and  the  Third  Mode  Backscatters 

21  or 

E  (t)  ■  T a  e  “  Cos(u>  t  +  $  ),  o<t<7.26. 


E  (t )  for  First-Mode  Excitation 


Natural  Mode  Components 


.2703  e~*5xCos(. 8661-91. 86°) 

2  .1085  e“'702TCos{l. 8071+128. 11°) 

3  .1245  e'‘843TCos(2. 7581+90. 38°) 

.1394  e  3  Cos(3.715r+46.48°) 

5  .1639  e"1-048TCos(4. 6761-5. 45°) 

-1  1 7Qr 

6  .2184  e  Cos(5.642t-62.5°) 


B 


.3265  e_1'20lTCos(6. 611-114. 95°) 
.4909  e-1‘267TCos(7. 581-161. 2°) 
.6914  e_1'327TCos (8. 551+155. 26°) 


10  .9013  e_1'382TCos(9. 5271+112. 82°) 

-1  434r 

11  1.108  e  Cos(10.5t+69.4°) 

12  1.304  e-1‘483TCos(ll. 481+22. 87°) 

—  1  *528r 

13  1.507  e  Cos(12. 461-26. 37°) 

14  1.839  e_1’572TCos(13. 441-78. 36°) 

15  2.337  e-1*613TCos(14. 421-128. 73°) 

16  3.010  e_1‘653TCos(15. 401-173. 66°) 

17  3.634  e-1*69lTCos(16. 381+148. 54°) 

-1  727t 

18  3.794  e  tCos(17. 36t+125.69°) 

—1  7A2r 

19  4.933  e  Xoa (18.341+156. 14°) 

20  9.203  e"2*°TCos(18. 8911-58. 56°) 

-  55t 

21  0.193  e  Cos(19. 501-79. 67°) 


E  for  Third-Mode  Excitation 


n  I  Natural  Mode  Components 


.2624  e”*5TCos(. 8661+171.1°) 

.2972  e"'702TCos (1.8071+138. 72°) 
2.743  e~'843TCos (2. 7581-83. 66°) 

_  Q  C£t 

.3245  e  *  Cos (3. 715T+68. 75°) 

—1  fi4ftr 

.2945  e  Cos(4.676x+21.09°) 

-1  129t 

.2880  e  Cos(5.642t-48.45q) 

.4556  e"1,2°lTCos(6. 611-120. 05°) 

.8078  e-1*267TCos(7. 581-172. 15°) 

1.264  e_1-327Tcos(8. 551+143. 81°) 
—1  3R2t 

10  1.752  e  Cos (9. 5271+102.09°) 

-1  434t 

11  2.233  e  Cos (10. 51+59. 56°) 

12  2.687  e_1'483TCos (11. 481+13. 75°) 

52Rt 

13  3.151  e  Cos(12. 46i-35. 07°) 

14  3.895  e_1*572TCos(13. 441-86. 90°) 

15  5.025  e“1’613TCos(14. 421-137. 11°) 

16  6.560  e_1‘653TCos(15. 401+178. 25°) 

17  8.011  e"1,691XCos(16. 381+140.79°) 

18  8.443  e"1,727TCos(17. 361+118. 21°) 

19  11.00  e'1,762TCos (18. 341+148. 72°) 

20  20.45  e-2,0xCos(18. 8911-65. 54°) 

_  ccT 

21  .0437  e  Co8(19.50t-86.75°) 


Relative  Amplitude 
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fig.  8  The  return  signals  from  the  right  sphere  and  a  wrong  sphere  which 
radius  is  10%  smaller  than  that  of  the  right  sphere  when  they  are 
illuminated  by  the  incident  signal  synthesized  to  excite  the  first 
natural  mode  of  the  right  sphere. 


fig.  9  The  return  signals  from  the  right  sphere  and  a  wrong  sphere  which 
radius  is  10X  smaller  than  that  of  the  right  sphere  when  they  are 
illuminated  by  the  incident  signal  synthesized  to  excite  the 
third  natural  mode  of  the  right  sphere. 
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1.  Work  in  Progress 

The  study  on  the  radar  waveform  synthesis  for  exciting  a  single¬ 
mode  backscatter  from  a  normally  oriented,  infinite  cylinder  is 
in  its  final  stage.  The  results  obtained  up  to  date  are  briefly 
outlined  below. 

A  radar  signal  propagating  in  the  4x-direction  is  incident  on 
an  perfectly  conducting,  infinite  cylinder  of  radius  a  with  its 
axis  located  along  the  z  axis.  The  electric  field  of  the  incident 
radar  signal,  in  its  Laplace  transform,  is  assumed  to  be 


Ei(r\ s) 


y  F<s)e^ae-Yrcos* 


(1) 


where  F(s)  is  an  unknown  waveform  function  to  be  synthesized  in 
such  a  way  that  If'Cr'.s)  excites  a  single-mode  backscatter  from 
the  cylinder,  and  y  ■  s/c. 

The  backscattered  electric  field  can  be  obtained,  through  a 
long  theoretical  development,  to  be 


Es  (r“,  s)  *  -y  F(s) 


*•.-*<*«  -  a>H(S) 


(2) 


where  r^  is  the  distance  between  the  observation  point  and  the 
cylinder  and  H(s)  is  the  transfer  function.  H(s)  can  be  expressed 
as 


H(s)  = 


-  l'n(0*ZK 

n=o  n  kJ,(C)VT 


(3) 


where  ln(C)  and  Kn(5)  are  the  derivatives  of  the  first  and  the 
second  kind  of  the  modified  Bessel  functions,  «  =1  for  n=0  and 
*2  for  n  >  0,  and  £ 


n 


ya  **  s(a/c)  . 

The  impulse  response  of  the  cylinder  is  obtained  by  inverting 
H(s) .  To  invert  H(s),  the  first  step  is  to  find  the  poles  of  H(s) 

•  t 

or  the  roots  of  K^C)  .  The  roots  of  Kn(£)  are  shown  in  Fig.  1. 

The  next  step  is  to  perform  the  integration  around  the  branch  cut 
of  H(s)  along  the  negative  real  axis  on  the  s-plane.  By  so  doing, 
the  impulse  response  of  the  cylinder  h(t)  can  be  expressed  as 
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e~e<t  '  2>In<»Kn<»  « 

■/ICK^U)  +  *2In2<£>] 

(4) 

t 

where  £  .  is  the  ith  root  of  the  K  (£)  function  and  a  .  is  the 
ni  n  ni 

residue  of  H(s)  at  Cn^-  The  first  term  of  eq.  (4)  is  the  sum  of 
all  the  natural  modes  and  the  second  is  a  line  integral  which  comes 
from  the  integration  around  the  branch  cut. 

Numerically ,  we  have  found  that  the  sum  of  a  finite  number 
of  natural  modes  plus  the  numerical  integration  of  the  line  inte¬ 
gral  of  eq.  (4)  given  an  approximate  impulse  response  which  is 
quite  accurate  for  the  late- time  period.  Preliminary  results  of 
the  approximate  impulse  response  of  the  cylinder  is  shown  in  Fig. 

2. 

We  are  in  the  process  of  calculating  the  required  incident 
signals  for  exciting  various  single-mode  backscatters .  Once  these 
required  incident  signals  are  obtained,  the  return  signals  can  be 
obtained  by  convoluting  the  incident  signals  with  the  impulse  res¬ 
ponse  of  the  cylinder.  We  expect  to  complete  these  calculations 
in  the  near  future. 


h(t) 


oo 

I 

n»o 


n  or 
n+1 

I 

i*l 


SniC 

Re(an.e  ni  ) 


oo 

z  if 

n=o  o 


2 .  Future  Plans 


The  following  topics  will  receive  major  attention  in  the  fu¬ 
ture  . 

1.  We  will  initiate  the  study  on  the  radar  waveform 
synthesis  for  exciting  single-mode  backscatters 
from  two  coupled  wire  targets  which  are  orientied 
at  an  angle  and  illuminated  by  an  incident  radar 
signal  at  an  arbitrary  direction.  The  case  of 
special  interest  will  be  a  wire  target  placed  over 
a  ground  plane  and  illuminated  by  an  incident  sig¬ 
nal.  Results  of  this  study  will  guide  the  experiment 
to  be  conducted  on  a  ground  plane . 


2.  We  will  initiate  the  study  on  the  cross-wire  structure 
which  simulates  an  aircraft. 

3.  We  will  design  an  experimental  setup  for  conducting 
the  experiment. 
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The  following  personnel  have  participated  in  this  research 
program. 

(1)  Kun-*fti  Chen,  Professor  and  principal  investigator. 

(2)  Dennis  P.  Nyquist,  Professor  and  senior  investigator. 

(3)  Byron  Drachman,  Associate  Professor  of  mathematics, 
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fig.  2. 


The  impulse  response  of  an  infinite  cylinder  of  radius  a 
when  the  impulse  signal  is  incident  normally  on  the  cylinder. 


